Staggered local density-of-states around the vortex in underdoped cuprates 
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We have studied a single vortex with the staggered flux (SF) core based on the SU(2) slave-boson 
theory of high Tc superconductors. We find that whereas the center in the vortex core is a SF 
state, as one moves away from the core center, a correlated staggered modulation of the hopping 
amplitude x s-nd pairing amplitude A becomes predominant. We predict that in this region, the 
local density-of-states (LDOS) exhibits staggered modulation when measured on the bonds, which 
may be directly detected by STM experiments. 



The elucidation of the ground state properties of the 
high Tc cuprates has been one of the major challenges 
in condensed matter physics. Recent STM experiments 
on vortices have indicated the existence of a normal 
core with a gap structure characteristic of the normal 
state pseudo gap above Tc- The theoretical description 
of the normal core, however, remains unresolved 1^-^. 
Based on the SU(2) slave-boson theory of the high Tc su- 
perconductors 1^,0, Lee and Wen ||^ proposed a model 
of the vortex with a staggered flux (SF) core, character- 
ized by a pseudo gap and staggered orbital current. The 
signature of the staggered flux (SF) order in zero field 
was found in the Gutzwiller-projected d-wave supercon- 
ducting state by variational Monte Carlo method ||] and 
exact diagonalization of the t-J model By using 

a Gutzwiller-projected U(l) slave-boson mean-field wave 
function, Han, Wang, and Lee recently found evidence of 
the SF order near the vortex core |11 1^. The natural 
question arises as to whether the SF core has observ- 
able consequences for the tunneling spectra using atomic 
resolution STM. This is the question addressed by the 
present paper. 

We first summarize some of the salient feature of the 
SU(2) vortex structure. In the SU(2) slave-boson the- 
ory |6|, the physical electron is represented by a spin- 
1/2 fermion operator fi^ and a charge-1 boson oper- 
ator hj = (feji,6i2): = -i=/i|Vio-, where hi and 

"i^Ia = Ui<y,Sasfls) are SU(2) doublets. The slave- 
boson mean-field state is characterized by 2 x 2 matri- 
ces C/„ - ^ "^'^ 
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on the links and boson con- 



densation hi on the sites. The Uij describe the hop- 
ping of the fermions and bosons. The mean-field so- 
lution is obtained by integrating out the fermions and 
minimizing the mean-field energy i?({?7y , /i^}), which 

leads to Xy = {fLfi<y) and A^^ = {easftafjs)- The 
SU(2) gauge invariance is realized through the relation 
E{{U,j,hi\) = E{{W,U,jWlWih,}), for any G 
SU(2). In the underdoped region, the mean-field solution 
has the following form (in the staggered flux (SF) gauge) 
t/f/ = -^r3exp[i(-l)*-+J«$oT-^] with A = + A2 
and <f>o — tan~^(A/x). This describes fermions hopping 



with flux ±4$o on alternating plaquettes [|T3|. The ad- 
vantage of the SF gauge is that it is apparent that the 
SU(2) symmetry has been broken down to the residual 
U(l), since Ufj' contains only . As a result, the low 
lying fluctuations have a simple form in the SF gauge: 

t/SF = -Ar^ exp [i(-l)'-+J»$oT^] exp [-ia%T^] . (1) 

These fluctuations are described by the lattice gauge field 

In the presence of a magnetic field, the mean-field so- 
lution contains vortices. Since the bosons are locally con- 
densed to the bottom of the band, they can be param- 
eterized by {h^^Y = ^/x{zii,—i{—\y'^^'^yzi2)^ where zn 
and Zi2 are slowly varying in space and time. The vortex 
is now described by the twisted Zi field and Uij in Eq. (|^). 
To write down the vortex structure explicitly, let us intro- 
duce Zii — e"^" cos Y and Zi2 = e*'^'^ sin where ipn = 
oti — 4>i/'^, and Lpi2 = ai+(l>i/2. The phase angles and (fii 
are associated with the electromagnetic (EM) U(l) and 
the U(l) gauge structures, respectively. The inter- 
nal degrees of freedom (f> and 9 can be visualized by the 
vector li — zjrzi = (sinff^ sin 0^, — sin^j cos cos 0^), 
which has the meaning of the quantization axis for the Zi 
fields. 

In the vortex structure proposed by Lee and Wen 
both ai and (pi/ 2 wind by tt and consequently give an 
appropriate hc/2e vortex for the EM gauge field A{r), 
i.e. Va = V(/)/2 = e^/2r which leads to Wipi — and 
V(/32 = e^/r, where denotes the azimuthal unit vector 
in the physical space. That is to say, only &2 changes its 
phase Lp2 by 27r as we go around the vortex, while bi does 
not. Outside the vortex core, li — (sin — cos 0^, 0) 
(see Fig. 1(a)), which describes a SC state. As we ap- 
proach the core, I62I must vanish and the vortex center 
is represented by li = (0, 0, 1), which is just the metaUic 
SF state. The /^-vector tilts smoothly from the equator 
to the north pole as the core is approached with a length 
scale denoted by ^c■ At the same time, the gauge flux 

= V X a'^ is distributed over a distance scale of Aa and 
it is expected that Aq ~ a;~^ > £c- Inside Aa, the gauge 
invariant combination 
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and its continuum limit, V(t") = iV(/)(r) — a^{r), is fi- 
nite. This is the analog of the superfluid velocity inside 
the London penetration depth Al in conventional vor- 
tices, except that here Aq ^ Al and the effect of the EM 
field A is negligible. 

The physics of the local electronic state is better visu- 
alized by making a local gauge transformation ipi^ — > 
Uff - g^U!f9l, and {hfV - {g^hfV = 
e'"' (-^/a;, 0), i.e. the local boson isospin vector points to- 
ward the north pole, as shown in Fig. 1(b). The explicit 

is given by gi = exp[i(— ^r^] exp[i^T^]. The 
advantage of this gauge is that the physical electron op- 
erator Cct c>c /cr. The new Uij has a physical meaning as 
governing the hopping and pairing of electrons. We shall 
refer to this as the c?-wave gauge. Indeed, locally there is 
a single boson component and the problem reduces to the 
more familiar U(l) mean-field theory, but with Xij and 
Aij which vary in space. This is precisely the problem 
treated by Han, Wang, and Lee [ pTp^ and it is grati- 
fying that they found numerically the staggered current 
around the vortex core as proposed in Ref. m. 




(a) (b) 
FIG. 1. (a) Configuration of the /i-vector outside the SF 
core in the SF gauge. At the center of the vortex, li points 
toward the north pole, which corresponds to the SF state. The 
local gauge transformation gi transforms this configuration to 
(b) in the d-wave gauge, where the internal phases of the hose 
condensate are gauged away. 

Now we are ready to ask the question: is any signature 
of the unit cell doubling directly measurable by STM 
tunneling? It turns out that there is no effect in the 
SF phase in the center of the core, because the period 
doubling of the current does not show up in the LDOS. 
This leads us to look outside the immediate core and 
examine the effect of the phase winding, i.e. the effect 
of Vij . After a local gauge transformation to the d-wave 
gauge, we find Eq. (Q) becomes 



Ufj = -X^3 \ COS 
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where 



and 



Xy = A cos , Ay = A sin , 



(4) 



(5) 



In this region, 9i ~ 9j ~ 7r/2 and Eq. (||) becomes 

Utj Xij'T^ + (-l)'''+-'"Ayr^ Recalling that Xij and 

Ay are interpreted as the hopping and pairing ampli- 
tudes of physical electrons, we see that the region outside 
the SF core is characterized by the staggered amplitude 
modulation. Note from Eq. that the amplitude of Xij 
and Aij are modulated in a correlated way to preserve 
Xij + Afj = constant. In Fig. 2, we schematically show 
the modulation pattern of Xij ■ In this region, Ufj breaks 
the translational symmetry, but does not break the time 
reversal symmetry. 




We now interpret this equation in different limits. 
First, we consider the region far outside the SF core. 



FIG. 2. Schematic drawing of the amplitude modulation 
pattern of the hopping field Xij outside the SF core. Solid 
and dotted bonds indicate enhanced and reduced amplitudes, 
respectively, where thickness of the bonds qualitatively repre- 
sents magnitude of the modulation. Circulation of the inter- 
nal gauge field, V(r), is indicated by the arrows. The dotted 
circle represents the SF core boundary, r — £c, inside which 
the staggered orbital current fiows. Note that the boundary 
should not be taken literally. In reality, there is a crossover 
region around r ~ where the staggered current and the 
staggered amplitude modulation coexist. 

Next, we consider the vicinity of the vortex cen- 
ter, where 0i ^ 9j ^ and Eq. (^) becomes Ufj ~ 
— ^r'^ exp l)*="+^«$y r'^] . In this region, Ufj breaks 
both the translational and the time reversal symmetries, 
and the staggered orbital current appears |^. How- 
ever, the staggered flux is further modulated according 
to Eq. 1). 

As we approach the core from the outside, the 7^ -vector 
in the SF gauge representation gradually rises off from 
the equatorial plane [see Fig. 1(a)]. This gives rise to a 
crossover region characterized by coexistence of the am- 
plitude and phase modulation, where the 0-dependence 
of Ufj becomes significant. However, to study the effects 
of a ^-dependent t/fj is beyond the scope of the present 
paper. From now on, we shall compute the LDOS by 
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setting 9i = 9j = n/2. We expect our results to be qual- 
itatively valid for r }^ £c as long as we avoid the exact 
core center. 

The presence of staggered modulation of Xij and A^^ 
in this region suggests that this may be the best place to 
look for unit cell doubling effects. To model the tunneling 
current, we assume that the electrons tunnel from the tip, 
located at r, to a linear combination of Wannier orbitals 
centered at lattice sites «, i.e. ~ '"O- Then 

the LDOS in the d-wave gauge is written as 

iV(r,w) = ^Im2_^ai(r)aj(r)[C/,?J-(iw)]ii \iuj=uj+iS ■ (6) 

hi 

The subscript 11 means the 11-component of the lattice 
fermion propagator in a 2 x 2 matrix form, Qfj (ioj) = 

Here we demonstrate that the LDOS exhibits conspicu- 
ous staggered pattern only when measured on the bonds. 
For example, we pick up the sites 1,2, ...,6 indicated in 
Fig. 2 and consider the midpoints on the bonds, Bi, B2, 
and the plaquette centers, Ci, C2. The LDOSs at Ci and 
C2 come from Eij=i,2,4,5 ^'^'^ Eij=2,3,5,6 ^ij ' respec- 
tively. We see, however, that ^56 since the bonds 
12 and 56 are almost equivalent except for a small in- 
equivalence coming from the nonuniformity of the slowly- 
varying V(r) field. Similarly, ^45 ^ Cjfg and Q^a, ^ Sfe- 
Therefore, N{Gi,uj) - iV(C2,w). Similarly, the LDOS 
at the lattice sites is almost uniform. On the other hand, 
the LDOSs at Bi and B2 come from -^-^ ^ Qf^ , and 

J2i j=2 3^fj' respectively. Here, G12 and f/fg are clearly 
inequivalent because they connect bonds with alternating 
hopping-pairing amplitudes. Furthermore, it is seen that 
the staggered modulation of the LDOS becomes most 
conspicuous when scanned along the a- or 6-axis [see the 
inset of Fig. 3(a)], because on these bonds the circulating 
V(r) field becomes parallel to the bond directions. From 
now on, we shall concentrate on the LDOS at the point 
r = {0,iy + 1/2) with lattice unit. To evaluate Qfj{iLLj), 
we shall use the following two approaches which may be 
complementary to each other: (I) gradient expansion, and 
(II) uniform V approximation. 
(I) Gradient expansion: 

First, we expand Eq. (^ with respect to Vij up to 
first order as Xij X— (— l)*"'^-'" AVy, and ~ A -|- 
(_l)i.+j„^Vy, which gives Ufj = Ufj+5Ufj with SUij = 
(_l)«x+i„[Ar3 + (-l)»«+i«;^ri]Vy. Then, we treat the 
effect of 5Uij within the Born approximation. In this 
case, we can take account of the circulating configuration 
of the V(r) field. 

The LDOS on the bonds is written as N{r,uj)/xa'^ 
— No{ll!) + SN{r,u!). The uniform counterpart is given 
by 7Vo(cj) = -^EfcCOS^^ Im[^J(fc, za;)]ii 
where gl{k,iuj) = Uk/iiuj - Ek) + Vk / {ioJ + Ek) , Uk = 
i[l + hkT^ + VkT^)/Ek], Vk^^[l- i-fkT^ + VkT^)/Ek], 



and Ek = + vl ^^^^ 7fe = -Jx[cosfca; -I- cosfcj^ -|- 
^2 cos fc^: cos fcy -I- t3{cos2kx + cos2fcj^)] -I- oq, and rjk = 
+JA{coskx — cos ky). We took account of the sec- 
ond and third nearest neighbor hopping of the fermions 
[^2 = —0.550 and is = 0.087] to reproduce the real band 
structure of Bi2Sr2CaCu208+5 as measured by angle- 
resolved photo emission spectroscopy Q . 

The perturbation term causes the unit cell doubling 
and is written in Fourier space as 

SH^ = -EfcGRzEq^smallE<x[V'I+^+Q^Cfc(q),/;fc_,^ + 

H.c.]. Recalling the expression for h{q), we obtain 
Ckiq) = AC+(q)r3 + xCkiq)r\ where C±(q) = 
TT J j^^^^j^ [qy sin kr^ ± sin ky] , RZ means the re- 
duced zone j/ca;! + |/cy| < tt, and Q — (tt, tt). We obtain 

SN{r,Lj) ^ t^Vxir) ^ sin" k, 
fceRZ 

[Ll6{uj; Ek, Ek+q) + L]^S{uj; -Ek, -Ek+q) 
-f N+5{uo; Ek,~Ek+q) + N^5{uo; -Ek, Ek+q)] , (7) 

where r = {0,iy + 1/2), 6{uj;x,y) = [S{lj — x) — S{iLj — 
y)]/{x-y), = A(l-l-7+7_-r/+r/_±7+±7_)±x(?7+ + 
11- ± 7+»7- ± ^ A(l - 7+7_ + r/+r/_ ± 7+ T 

J-)±X{V+-V-T1+V-TV+J-), 7+ = Ik+q/Ek+q, 7- = 
Ik/Ek, r]+ = rjk+q/Ek+q, and 77. =r]k/Ek. Because of 
the staggerdness, a fermion with k in the reduced zone 
is scattered to fc -I- Q in the second zone. 

In Fig. 3(a), we show the profile of N {r , uj) / xa^ at the 
points A{0, iy + 1/2), B(l/2,iy + 1), C{0, iy + 3/2), and 
D(— 1/2, iy + 1) by choosing iy to satisfy Vx{r) — 0.1. As 
a rough estimate, this corresponds to a lattice position 
of iy — A if Xa — 10. We used oq = 0.05%^ and A/x — 
0.2. Numerical integration was performed by dividing the 
Brillouin zone into a 320 x 320 mesh. Note that at B and 
D, SN{r,uj) almost vanishes and the LDOS is just given 
by Nq(uj), because V(t) becomes almost perpendicular 
to these bond directions. The modulation pattern at the 
other points can be read off from Fig. 2. 

It is remarkable that, inside the overall V-shaped pro- 
file with the sharp peaks at d) = oj/xJ = ±0.323 asso- 
ciated with the c?-wave superconducting gap, there ap- 
pears additional staggered structure around Co = ±0.179 
and uj = 0.226. This structure comes from resonant scat- 
tering between the fermions with k and k + Q. As lo 
increases from zero, the energy contours Ek = w and 
Ek+q = LO first touch on the reduced zone boundary 
at ci) = ±0.179, as indicated in Fig. 3(b), and reso- 
nance occurs. Then, at co — ±0.226, they touch again at 
{tt/2, 7r/2) and the second resonance occurs. The second 
resonance comes up only in the electron {to > 0) side due 
to the matrix element effect [L^ vanishes at (7r/2, 7r/2)]. 
It is naturally expected that this resonant scattering may 
open a gap in the fermion excitation spectrum if we go 
beyond the perturbative scheme. This point can be con- 
firmed through the exact treatment under the uniform V 
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approximation shown below. 
(II) Uniform V approximation: 

Next we consider the case of uniform V = {Vx,Vy) 
which may locaUy capture the effects of the circulating 
V(r). In this case, we exactly diagonalize the fermion 
Hamiltonian. The LDOS on the bonds is written in a 
form N{uj)/ xo? = iVo(w) ± ^/^(a;), where ± signs al- 
ternate from bond to bond. As was inferred from the 
perturbative analysis, the unit cell doubling splits the 
one-particle spectrum into two branches [see Fig. 3(e)], 
±El and ±E]^. 




FIG. 3. (a) LDOS profile at the points A, B, C, and D indi- 
cated in the inset. The LDOS at B and C are just No{ui). The 
energy contours Ek = uj and Ek+Q = u) touch at (5 = 0.179 
and UJ = 0.226, as indicated in (b) and (c), respectively. The 
sharp peaks at = ±0.323 are associated with the supercon- 
ducting gap. The small wiggles outside the V-shaped profile 
come from numerical fluctuations, (d) Profile of No{ui), and 
No{uj) ± 5N{uj) for the uniform gauge fleld V = (-0.1, 0). (e) 
The band dispersion of the splitted bands iJ^ along the path 
r(0,0) M(7r/2,7r/2) Y(0, tt) ^ F. Fine band splittings 
on the reduced zone boundary are magnifled in the inset, (f) 
Fine structure of No{u>) around ui ~ 0.2, detected with higher 
numerical resolution. 



In Fig. 3(d), we show the profile of No{^) and iVo(a;)± 
SN{uj) for V = (—0.1,0), the direction and strength of 
which correspond to V(t) around the points B, D, and 
A,C in Fig. 3(a), respectively. We used the same pa- 
rameter set as in the case of Fig. 3(a). We see that the 
staggerd modulation profile, No{uj)±SN{uj), is in remark- 
able agreement with that obtained by the perturbative 
analysis. The structures correspond to van-Hove singu- 
larities associated with the gap opening on the reduced 
zone boundary as shown in Fig. 3(e). However, a striking 
difference is that the dip structure around uj = 0.2 ap- 
pears even in the uniform counterpart iVo(a-'). This sug- 
gests that in reality the dip structure may be detected 



not only on the bonds but also at sites. In Fig. 3(d), 
due to numerical resolution [320x320 mesh of the Bril- 
louin zone] we resolve structures around lu — ±0.2 and 
UJ = 0.26. As indicated in Fig. 3(e), there are fine band 
splittings on the reduced zone boundary leading to van 
Hove singularities at <I> = 0.186, 0.216, 0.229, and 0.265. 
The corresponding fine structure in Nq(lu) could be de- 
tected with much higher numerical resolution [720x720 
mesh of the Brillouin zone], as shown in Fig. 3.(f). 

We note that, in both approximations, the modulated 
structure in the LDOS is predominant on the particle 
side {oj > 0). We can understand this asymmetry by 
first turning off the superconductivity and considering 
the effect of unit cell doubling. Since we are doping with 
holes, the gaps being opened by the unit cell doubling are 
on the empty side of the Fermi surface. Matrix element 
effects preserve this particle-hole asymmetry even after 
we turn on the superconductivity. 

Combining the results obtained through the gradient 
expansion and the uniform V approximation, we may 
reasonably say that the signature of the unit cell doubling 
may be most prominently detected through the charac- 
teristic dip structure inside the V-shaped profile. The 
structure predicted here is very specific and its observa- 
tion will be a strong confirmation of the SU(2) vortex 
model. 
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